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Gravitational torques between a planet and gas in the protoplanetary disk result in orbital mi- 
gration of the planet, and modification of the disk surface density. Migration via this mechanism 
is likely to play an important role in the formation and early evolution of planetary systems. 
For masses comparable to those of observed giant extrasolar planets, the interaction with the 
disk is strong enough to form a gap, leading to coupled evolution of the planet and disk on a 
viscous time scale (Type II migration). Both the existence of hot Jupiters, and the statistical 
distribution of observed orbital radii, are consistent with an important role for Type II migration 
in the history of currently observed systems. We discuss the possibility of improving constraints 
on migration by including information on the host stars' metallicity, and note that migration 
could also form a population of massive planets at large orbital radii that may be indirectly 
detected via their influence on debris disks. For lower mass planets with Mp ~ M^, surface den- 
sity perturbations created by the planet are small, and migration in a laminar disk is driven by 
an intrinsic and apparently robust asymmetry between interior and exterior torques. Analytic 
and numerical calculations of this Type I migration are in reasonable accord, and predict rapid 
orbital decay during the final stages of the formation of giant planet cores. The difficulty of 
reconciling Type I migration with giant planet formation may signal basic errors in our under- 
standing of protoplanetary disks, core accretion, or both. We discuss physical effects that might 
alter Type I behavior, in particular the possibility that for sufficiently low masses (Mp — > 0) 
turbulent fluctuations in the gas surface density dominate the torque, leading to random walk 
migration of very low mass bodies. 



1. Introduction 

The extremely short orbital period of 51 Pegasi (Mayor & Queloz 1995) and the other 
hot Jupiters pose a problem for planet formation, not only because such systems bear 
little resemblance to the Solar System, but more fundamentally because the high tem- 
peratures expected in the protoplanetary disk at radii a < 0.1 AU largely preclude the 
possibility of in situ formation. Disk models by Bell et al. (1997) show that for typical 
T Tauri accretion rates of M ~ 10~* MQyr~^ (GuUbring et al. 1998) the midplane tem- 
perature interior to 0.1 AU exceeds 1000 K, destroying ices and, for the very closest in 
planets, even dust. At least the cores of these hot Jupiters must therefore have formed 
elsewhere, and subsequently migrated inward. Migration is also likely to have occurred 
for the larger population of extrasolar planets that now lie within the snow line in their 
parent disks (Bodenheimer, Hubickyj & Lissauer 2000), though this is a more model- 
dependent statement since both the location of the snow line (Sasselov & Lecar 2000) 
and its significance for giant planet formation remain uncertain. 

Orbital migration of planets involves a loss of angular momentum to either gas or 
other solid bodies in the system. Three main mechanisms have been proposed, all of 
which involve purely gravitational interactions (aerodynamic drag, which is central to 
the orbital evolution of meter-scale rocks, is negligible for planetary masses). The first 
is gravitational interaction between the planet and the gas in the protoplanetary disk. 
This leads to angular momentum exchange between the planet and the gas, and resulting 
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orbital evolution (Goldreich & Tremaine 1980; Lin, Bodenheimer & Richardson 1996). 
Since gas giant planets, by definition, formed at an epoch when the protoplanetary disk 
was still gas-rich, this type of migration is almost unavoidable. It is the main subject of 
this article. However, further migration could also occur later on, after the gas disk has 
been dissipated, as a consequence of the gravitational scattering of either planetesimals 
(Murray et al. 1998) or other massive planets (Rasio & Ford 1996; Weidenschilling & 
Marzari 1996; Lin & Ida 1997; Papaloizou & Terquem 2001). Some orbital evolution 
from planetesimal scattering is inevitable, given that the formation of massive planets is 
highly likely to leave a significant mass of smaller bodies in orbits close enough to feel 
perturbations from the newly formed giant. In the Solar System, planetesimal scattering 
could have allowed substantial outward migration (Thommes, Duncan & Levison 1999) 
of Uranus and Neptune - which have a small fraction of the Solar System's angular 
momentum - while simultaneously raising the eccentricities and inclinations of all the 
giant planets to values consistent with those observed (Tsiganis et al. 2005). Although 
this is an attractive theory for the architecture of the outer Solar System, invoking a 
scaled-up version of this process as the origin of the hot Jupiters is problematic. To drive 
large-scale migration of the typically rather massive planets seen in extrasolar planetary 
systems would rciquirc a comparable mass of planetesimals interior to the initial orbit of 
the planet. Such a planetesimal disk would in turn imply the prior existence of a rather 
massive gas disk, which would likely be more effective at causing migration than the 
planetesimals. Similar reservations apply to models of planet-planet scattering, which is 
only able to yield a population of planets at small orbital radii if multiple planet formation 
(with the planets close enough that they are unstable over long periods) is common. That 
said, the observation that most extrasolar planets have significantly eccentric orbits - 
which currently defies explanation except as an outcome of planet-planet scattering (Ford, 
Rasio & Yu 2003) - may mean that at least some scattering-driven migration occurs in 
the typical system. 

Figure 1 illustrates how a planet on a circular orbit interacts with the protoplanetary 
disk. The planet perturbs the gas as it passes by the planet, with angular momentum 
transport taking place at the locations of resonances in the disk - radii where a char- 
acteristic disk frequency is related to the planet's orbital frequency. For relatively low 
mass perturbers, the interaction launches a trailing spiral wave in the gas disk, but is not 
strong enough to significantly perturb the azimuthally averaged surface density profile. 
In this regime, described as Type I migration, angular momentum transport between the 
planet and the gas occurs while the planet remains embedded within the protoplanetary 
disk. The rate of migration is controlled by the sum of the torques arising from the inner 
and outer Lindblad and corotation resonances, which is generally non-zero (if the sum 
happened to be close to zero, the planet would act as a source of angular momentum 
transport in the disk (Goodman & Rafikov 2001) while remaining in place). For the 
parameters (sound speed, efficiency of angular momentum transport) that are believed 
to be appropriate for protoplanetary disks. Type I migration occurs for planet masses 
Mp < 0.1 Mj, where Mj is the mass of Jupiter (Bate et al. 2003), and is most rapid 
as this critical mass is approached (e.g. Ward 1997). As a result, it is likely to play a 
particularly important role in the final assembly of giant planet cores. 

At higher masses - Mp > 0.1 Mj the angular momentum removal / deposition 
at the planet's inner / outer Lindblad resonances is strong enough to repel gas from 
an annular region surrounding the planet's orbit, forming a gap in which the surface 
density is reduced compared to its unperturbed value. For planets of a Jupiter mass 
and above, the gap is almost entirely evacuated (e.g. the right-hand panel of Figure 1), 
although mass may continue to flow in a stream across the gap to enter the planet's Roche 
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Figure 1 . An illustration of the interaction between a planet on a fixed circular orbit with a 
laminar (non-turbulent) protoplanetary disk, computed from a two-dimensional (r, <f)) hydrody- 
namic simulation with a locally isothermal equation of state and a constant kinematic viscosity. 
In the left-hand panel showing the regime of Type I migration, a relatively low mass planet 
excites a noticeable wave in the disk gas but does not significantly perturb the azimuthally aver- 
aged surface density profile (shown as the inset graph). In contrast, a 10 Mj planet (right-hand 
panel) clears an annular gap in the disk, within which the surface density is a small fraction 
of its unperturbed value. As the disk evolves over a viscous time scale, the planet is predicted 
to track the motion of the gas (either inward or outward) while remaining within the gap. 
This is Type II migration. For a movie showing the transition between these regimes, go to 
http : //jilawww. Colorado . edu/ ~pja/planet_migration .html. 



lobe (Artymowicz & Lubow 1996; Lubow, Seibert & Artymowicz 1999). The location of 
the inner and outer edges of the gap are set by a balance between angular momentum 
exchange with the planet (which tends to widen the gap), and internal stresses within 
the protoplanetary disk ('viscosity', which tends to close it). Under most circumstances, 
this balance acts to lock the planet into the long term viscous evolution of the disk gas. 
At radii where the disk gas is moving inward the planet migrates toward the star, all the 
while remaining within its gap (Lin & Papaloizou 1986). This is Type II migration, which 
differs from Type I not only in the presence of a gap, but also because the rate depends 
directly on the efficiency of angular momentum transport within the protoplanetary disk. 

In addition to these well-established migration regimes, qualitatively different behav- 
ior may occur at very low masses, and at masses intermediate between the Type I and 
Type II regimes. At sufficiently low masses (probably of ^ 10 and below) the per- 
sistent Type I torque may be overwhelmed by random torques from surface density 
perturbations in a turbulent protoplanetary disk (Nelson & Papaloizou 2004; Laughlin, 
Steinacker & Adams 2004). This process, which is similar to the heating of Galactic stars 
by transient spiral arms (Carlberg & Sellwood 1985), may lead to rapid random walk 
migration on top of Type I drift. This could have important consequences for both core 
accretion (Rice & Armitage 2003) and terrestrial planet formation. Another uncertain 
regime lies at the transition mass {Mp w 0.1 Mj) between Type I and Type II migra- 
tion, where a partial gap exists and corotation torques can be highly significant. Masset 
& Papaloizou (2003) suggested that the corotation torques could drive an instability in 
the direction of migration, which if confirmed would be extremely important for massive 
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planet formation. Subsequent higher resolution simulations by D'Angelo, Bate & Lubow 
(2005) have demonstrated that extremely high resolution is needed in order to accurately 
capture torques arising from within and around the Hill sphere. For a Saturn mass planet, 
the resolution had to be increased to approximately 50 grid zones per Hill radius before 
a numerically converged solution for the migration rate was attained. Lower resolutions 
yielded artificially rapid migration. D'Angelo, Bate & Lubow (2005) conclude from this 
that although torques from the near-planet region are significant, they do not appear to 
dramatically increase the migration rate. 



2. Type I migration 

In the Type I regime the perturbation induced by the planet in the gas disk remains 
small, and the net torque has surprisingly little dependence on the microphysics of the 
protoplanctary disk (e.g. Goldreich & Tremaine 1978; Lin & Papaloizou 1979). In par- 
ticular, viscosity - normally the most uncertain element of a protoplanctary disk model 
- enters only indirectly via its influence on the magnitude and radial gradient of the 
surface density and sound speed. Gencrically, the net torque scales with the planet mass 
as T oc Mp, so that the migration time scale at a given radius scales as r oc M~^. Type I 
migration therefore becomes increasingly important as the planet mass increases, and is 
fastest just prior to gap opening (the onset of which does depend on the disk viscosity). 
Despite its attractive lack of dependence on uncertain disk physics, the actual calcula- 
tion of the net torque is technically demanding, and substantial improvements have been 
made only recently. Here, we summarize a few key results - the reader is directed to the 
original papers (primarily by Artymowicz, Ward, and their collaborators) for full details 
of the calculations. 

2.1. Analytic calculations 
The simplest calculation of the Type I torque (Goldreich & Tremaine 1978, 1979, 1980) 
neglects significant pressure effects in the disk close to the planet, and is therefore valid 
for low m resonances. In this approximation, Lindblad resonances occur at radii in the 
disk where the epicyclic frequency n is an integral multiple m of the angular velocity 
in a frame rotating with the planet at angular velocity ilp. For a Keplerian disk, this 
condition, 

D{r) = K^-m'^{VL- flpf = 0, 
can be simplified using the fact that k = fl. The resonances lie at radii. 



rL={l±-) 
m I 



2/3 



where Vp is the planet's orbital radius. The lowest order resonances lie at r = 1.587rp 
and r = 0.630rp, but an increasingly dense array of high m resonances lie closer to the 
planet. Resonances at r < add angular momentum to the planet, while those at r > 
remove angular momentum. The torque at each resonance can be evaluated in terms 
of a forcing function as. 



rdD/dr 

where E is the gas surface density. Explicit expressions for are given by Goldreich & 
Tremaine (1979). The net torque is obtained by evaluating the torque at each resonance, 
and then summing over all m. 
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For Type I migration the behavior of gas close to the planet, where is largest, is 
critical. Accurately treating this regime requires elaborations of the basic Goldreich & 
Tremaine (1979) approach. For a razor thin two-dimensional disk model, - the effects of 
radial pressure and density gradients, the calculation is described in Ward (1997), and 
references therein (especially Ward 1988; Artymowicz 1993a, 199b; Korycansky & Pollack 
1993). These papers include the shifts in the location of Lindblad resonances due to both 
radial and azimuthal pressure gradients, which become significant effects at high m. For 
Lindblad resonances, the result is that the dominant torque arises from wavenumbers 
m ^ Tp/h^ where /i, the vertical disk scale height, is given in terms of the local sound 
speed Cs by h = Cs/ftp. The fractional net torque 2|rinner + router|/(|7inncr| + |routcr|) 
can be as largo as 50% (Ward 1997), with the outer resonances dominating and driving 
rapid inward migration. Moreover, the small shifts in the locations of resonances that 
occur in disks with different radial surface density profiles conspire so that the net torque 
is only weakly dependent on the surface density profile. This means that the predicted 
rapid inward migration occurs for essentially any disk model in which the sound speed 
decreases with increasing radius (Ward 1997). Corotation torques - which vanish in the 
oft-considered disk models with S oc r~^/^ - can alter the magnitude of the torque but 
are not sufficient to reverse the sign (Korycansky & Pollack 1993; Ward 1997). 

The observation that the dominant contribution to the total torque comes from gas 
that is only Ar ~ h away from the planet immediately implies that a two-dimensional 
representation of the disk is inadequate, even for protoplanetary disks which are geomet- 
rically thin by the usual definition {h/r < 0.1, so that pressure gradients, which scale as 
order (h/r)"^, are sub-percent level effects). Several new physical effects come into play 
in a three-dimensional disk: 

1. The perturbing potential has to be averaged over the vertical thickness of the disk, 
effectively reducing its strength for high m resonances (Miyoshi et al. 1999). 

2. The variation of the scale height with radius decouples the radial profile of the 
midplane density from that of the surface density. 

3. Wave propagation in three-dimensional disks is fundamentally different from that 
in two dimensions, if the vertical structure of the disk departs from isothermality (Lubow 
& Ogilvie 1998; Bate et al. 2002). 

Tanaka, Takeuchi & Ward (2002) have computed the interaction between a planet and 
a three-dimensional isothermal disk, including the first two of the above effects. Both 
Lindblad and corotation torques were evaluated. They find that the net torque is reduced 
by a factor of 2-3 as compared to a corresponding two-dimensional model, but that 
migration remains inward and is typically rapid. Specifically, defining a local migration 
time scale via, 

-rp 

Tanaka, Takeuchi & Ward (2002) find that for a disk in which S oc r~^, 



As expected, the time scale is inversely proportional to the planet mass and the local 
surface density. Since the bracket is ~ (/i/r^)^, the time scale also decreases quite rapidly 
for thinner disks, reflecting the fact that the peak torque arises from closer to the planet 
as the sound speed drops. 

Although still limited by the assumption of isothermality, the above expression repre- 
sents the current 'standard' estimate of the Type I migration rate. Although slower than 
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two-dimensional estimates, it is still rapid enough to pose a potential problem for planet 
formation via core accretion. There is therefore interest in studying additional physical 
effects that might reduce the rate further. The influence of disk turbulence is discussed 
more fully in the next Section; here we list some of the other effects that might play a 
role: 

1. Realistic disk structure models. The run of density and temperature in the 
midplane of the protoplanetary disk is not a smooth power-law, due to sharp changes in 

opacity and, potentially, the efficiency of angular momentum transport (Gammic 1996). 
Menou & Goodman (2004) have calculated Type I rates in Shakura-Sunyaev type disk 
models, and find that even using the standard Lindblad torque formula there exist regions 
of the disk where the migration rate is locally slow. Such zones could be preferred sites 
of planet formation. 

2. Thermal effects. Jang-Condell & Sasselov (2005) find that the dominant non- 

axisymmctric thermal effect arises from changes to the stellar illumination of the disk 
surface in the vicinity of the planet. This effect is most important at large radii, and can 
increase the migration time scale by up to a factor of 2 at distances of a few AU. 

3. Wave refiection. The standard analysis assumes that waves propagate away from 
the planet, and are dissipated before they reach boundaries or discontinuities in disk 
properties that might reflect them back toward the planet. Tanaka, Takeuchi & Ward 
(2002) observe that reflection off boimdarics has the potential to substantially reduce the 
migration rate. We note, however, that relaxation of vertical isothermality will probably 
lead to wave dissipation in the disk atmosphere within a limited radial distance (Lubow 
& Ogilvic 1998), and thereby reduce the possibility for reflection. 

4. Accretion. Growth of a planet during Type I migration is accompanied by a non- 
resonant torque, which has been evaluated by Nelson & Benz (2003a). If mass is able to 
accrete freely on to the planet, Bate et al. (2003) find from three-dimensional simulations 
that Mp (X Mp for Mp < 10 Mq, with a mass doubling time that is extremely short (less 
than 10^ yr). In reality, it seems likely the planet will be unable to accept mass at such 
a rapid rate, so the mass accretion rate and resulting torque will then depend on the 
planet structure. 

5. Magnetic fields. The dominant field component in magnetohydrodynamic disk 

turbulence initiated by the magnctorotational instability is toroidal (Balbus & Hawlcy 
1998). Terquem (2003) finds that gradients in plausible toroidal magnetic fields can signif- 
icantly alter the Type I rate, and in some circumstances (when the field decreases rapidly 
with r) stop migration. More generally, a patchy and variable toroidal field might lead 
to rapid variations in the migration rate. Whether this, or density fiuctuations induced 
by turbulence, is the primary influence of disk fields is unclear. 

6. Multiple planets. The interaction between multiple planets has not been studied 
in detail. Thommes (2005) notes that low mass planets, which would ordinarily suffer 
rapid Type I migration, can become captured into resonance with more massive bodies 
that arc themselves stabilized against rapid decay by a gap. This may be important 
for understanding multiple planet formation (and we have already noted that there is 
circumstantial evidence that multiple massive planet formation may be common), though 
it does not explain how the first planet to form can avoid rapid Type I inspiral. 

7. Disk Eccentricity. Type I migration in an axisymmetric disk is likely to damp 
planetary eccentricity. However, it remains possible that the protoplanetary disk itself 
might be spontaneously unstable to development of eccentricity (Ogilvic 2001). Pa- 
paloizou (2002) has shown that Type I migration can be qualitatively altered, and even 
reversed, if the background flow is eccentric. 
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2.2. Numerical simulations 

Hydrodynamic simulations of the Type I regime within a shearing sheet geometry have 
been presented by Miyoshi et al. (1999), and in cyHndrical geometry by D'Angelo, Hen- 
ning & Klcy (2002), D'Angelo, Klcy & Hcnning (2003), and Nelson & Bcnz (2003b), 
with the latter paper focusing on the transition between Type I and Type II behavior. 
The most comprehensive work to date is probably that of Bate et al. (2003), who sim- 
ulated in three dimensions the interaction with the disk of planets in the mass range 
1 M© ^ ^ 1 Mj. The disk model had h/r = 0.05, a Shakura-Sunyaev (1973) a 
parameter a = 4 x 10~^, and a fixed locally isothermal equation of state (i.e. Cg = Cs{r) 
only). This setup is closely comparable to that assumed in the calculations of Tanaka, 
Takeuchi & Ward (2002), and very good agreement was obtained between the simula- 
tion results and the analytic migration time scale. Based on this, it seems reasonable to 
conclude that within the knoiun limitations imposed by the restricted range of included 
physics^ current calculations of the Type I rate are technically reliable. Given this, it is 
interesting to explore the consequences of rapid Type I migration for planet formation 
itself. 

2.3. Consequences for planet formation 

The inverse scaling of the Type I migration time scale with planet mass means that the 
most dramatic effects for planet formation occur during the growth of giant plant cores via 
core accretion (Mizuno 1980). In the baseline calculation of Pollack et al. (1996), which 
does not incorporate migration, runaway accretion of Jupiter's envelope is catalyzed 
by the slow formation of a ~ 20 core over a period of almost 10 Myr. This is in 
conflict with models by Guillot (2004), which show that although Saturn has a core of 
around 15 M©, Jupiter's core is observationally limited to at most « 10 Mq, leading 
to discussion at this meeting of several ways to reduce the theoretically predicted core 
mass. Irrespective of the uncertainties, howcivcr, it seems inevitable that planets forming 
via core accretion pass through a relatively slow stage in which the growing planet has a 
mass of 5 — 10 M©. This stage is vulnerable to Type I drift. 

Figure 2 shows the migration time scale for a 10 M© planet within gas disks with 
surface density profiles of E oc (very roughly that suggested by theoretical models, 
e.g. Bell et al. 1997) and S oc r~^/^ (the minimum mass Solar Nebula profile of Weiden- 
schilling 1977). We consider disks with integrated gas masses (out to 30 AU) of 0.01 Mq 
and 1 Mj. The latter evidently represents the absolute minimum gas mass required to 
build Jupiter or a typical extrasolar giant planet. The torque formula of Tanaka, Takeuchi 
& Ward (2002) is used to calculate the migration time scale t. As is obvious from the 
Figure, migration from 5 AU on a time scale of 1 Myr - significantly less than either the 
typical disk lifetime (Haisch, Lada & Lada 2001) or the duration of the slow phase of 
core accretion - is inevitable for a core of mass 10 M^, even if there is only a trace of gas 
remaining at the time when the envelope is accreted. For more reasonable gas masses, the 
typical migration time scale at radii of a few AU is of the order of 10^ yr. Another repre- 
sentation of this is to note that in the giant planet forming region, wc predict significant 
migration (t = 10 Myr) for masses Mp > 0.1 M©, and rapid migration (r = 1 Myr) for 
Mp > 1 Mg. We can also plot, for the same disk model, an estimate of the isolation mass 
(Lissauer 1993; using the gas to planctesimal surface density scaling of Ida & Lin 2004a). 
The isolation mass represents the mass a growing planet can attain by consuming only 
those planetesimals within its feeding zone - as such it is reached relatively rapidly in 
planet formation models. Outside the snow line, the migration time scale for a planet 
(or growing core) at the isolation mass is typically a few Myr, reinforcing the conclusion 
that migration is inevitable in the early stage of giant planet formation. By contrast, in 
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Figure 2. Left hand panel: planet mass for which the Type I migration time scale at different 

radii equals 1 Myr (upper solid curve) and 10 Myr (lower curve). The gas disk is assumed to have 
a surface density profile E oc r^'^^^, and a mass within 30 AU of 0.01 Mq. The dotted lino shows 
an estimate of the isolation mass in the same disk model, assuming Solar motallicity and a snow 
lino at 2.7 AU. Right hand panel: the migration time scale for a 10 core in protoplanetary 
disks with surface density profiles of S oc r~^^^ (solid lines) and S oc (dashed lines). For 
each model, the lower curve shows the time scale for a disk with a gas mass of 0.01 M© within 
30 AU, while the upper curve shows results for an absolute minimum mass gas disk with only 
1 Mj within 30 AU. 



the terrestrial planet region, interior to the snow line, planets need to grow significantly 
beyond isolation before rapid migration ensues. It is therefore possible for the early stages 
of terrestrial planet formation to occur in the presence of gas with only limited Type I 
migration, while the final assembly of terrestrial planets happens subsequently in a gas 
poor environment. 

Does Type I migration help or hinder the growth of giant planets via core accretion? 

This question remains open, despite a history of investigations stretching back at least 
as far as papers by Hourigan & Ward (1984) and Ward (1989). Two competing effects 
are at work: 

1 . A migrating core continually moves into planetesimal-rich regions of the disk that 
have not been depleted by the core's prior growth. This depletion is, in part, responsible 
for the slow growth of Jupiter in static core calculations. Calculations suggest that a 
rapidly migrating core can capture of the order of 10% of the planetesimals it encoim- 
ters (Tanaka &: Ida 1999), with the collision fraction increasing with migration velocity. 
Slow migration velocities allow for inward shepherding of the planetesimals rather than 
capture, and a low accretion rate (Ward & Hahn 1995; Tanaka & Ida 1999). 

2. A migrating core must reach the critical core mass before it is lost to the star, on 
a time scale that, as we noted above, can be an order of magnitude or more smaller 
than the gas disk lifetime. Unfortunately, a high accretion rate of planetesimals increases 
the critical core mass needed before runaway gas accretion starts and, at fixed accretion 
rate, the critical mass also increases as the core moves inward (Papaloizou & Terquem 
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1999). Migration therefore favors a high rate of planetesimal accretion, but often hinders 
attaining the critical core mass needed for envelope accretion. 

Calculations of giant planet formation including steady core migration have been pre- 
sented by several groups, including Papaloizou & Terqucm (1999), Papaloizou & Larwood 
(2000), AUbcrt, Mordasini & Benz (2004) and Alibcrt ot al. (2005). The results suggest 
that, for a single growing core, Type 1 migration at the standard rate of Tanaka, Takeuchi 
& Ward (2002) is simply too fast to allow giant planet formation to occur across a reason- 
able range of radii in the protoplanetary disk. Most cores are lost to the star or, if they 
manage to accrete envelopes at all, do so at such small radii that their ultimate survival 
is doubtful. More leisurely migration, on the other hand, at a rate suppressed from the 
Tanaka, Takeuchi & Ward (2002) value by a factor of 10 to 100, helps core accretion by 
mitigating the depletion effect that acts as a bottleneck for a static core (Alibert et al. 
2005). 

This difficulty in reconciling our best estimates of the Type 1 migration rate with 
core accretion signals that something is probably wrong with one or both of these the- 
ories. Three possibilities suggest themselves. First, the Type I migration rate may be 
a substantial overestimate, by an order of magnitude or more. If so, there is no need 
for substantial changes to core accretion theory or to protoplanetary disk models. We 
have already enumerated a list of candidate physical reasons for why the Type I rate 
may be wrong, though achieving a sufficiently large suppression does not seem to be 
straightforward. Second, angular momentum transport may be strongly suppressed in 
the giant planet formation region by the low ionization fraction, which suppresses MHD 
instabilities that rely on coupling between the gas and the magnetic field (Gammie 1996; 
Sano et al. 2000). An almost inviscid disk could lower the threshold for gap opening 
sufiiciently far that the slow stage of core accretion occurred in a gas-poor environment 
(elements of such a model have been explored by Rafikov 2002; Matsumara & Pudritz 
2005). Finally, and perhaps most attractively, it may be possible to find a variant of core 
accretion that is compatible with undiluted Type 1 migration. For a single core, we have 
studied simple models in which random walk migration leads to large fluctuations in the 
planetesimal accretion rate and an early onset of criticality (Rice & Armitage 2003). In 
the more realistic situation where multiple cores are present in the disk, it is possible 
that the early loss of the first cores to form (at small radii just outside the snow line) 
could evacuate the inner disk of planetesimals, allowing subsequent cores to reach their 
critical mass and accrete envelopes as they migrate inward. Further work is needed to 
explore such scenarios quantitatively. 



3. Stochastic migration in turbulent disks 

To a first approximation the efficiency of angular momentum transport (unless it is very 
low) has little impact on the predicted Type I migration rate. This assumes, however, that 
the disk is laminar. More r(^alistically, angular momentum transport itself derives from 
turbulence, which is accompanied by a spatially and temporally varying pattern of density 
fluctuations in the protoplanetary disk. These fluctuations will exert random torques on 
planets of any mass embedded within the disk, in much the same way as transient spiral 
features in the Galactic disk act to increase the velocity dispersion of stellar populations 
(Carlberg & Sellwood 1985). If we assume that the random torques are uncorrelated with 
the presence of a planet, then the random torques' linear scaling with planet mass will 
dominate over the usual Type I torque (scaling as M^) for sufficiently low masses. The 
turbulence will then act to increase the velocity dispersion of coUisionless bodies, or, in 
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Figure 3. Structure in a turbulent disk accretion flow, from a global cylindrical MHD simula- 
tion (after Armitage 1998). The flow is here visualized using density (left hand panel) and the 
vertically averaged Maxwell stress (BrB^) (right hand panel) as tracers. In magnetically active 
disks, transient density and magnetic field fluctuations are present across a wide range of spatial 
scales. 

the presence of damping, to drive a random walk in the semi-major axis of low mass 
planets. 

To go beyond such generalities, and in particular to estimate the crossover mass be- 
tween stochastic and Type I migration, we need to specify the source of turbulence in 
the protoplanetary disk. MHD disk turbulence (see Figure 3 for an illustration) driven 
by the magnetorotational instability (Balbus & Hawley 1998) provides a well-understood 
source of outward angular momentum transport in sufficiently well-ionized disks, and has 
been used as a model system for studying stochastic migration by Nelson & Papaloizou 
(2004) and by Laughlin, Steinacker & Adams (2004). Density fluctuations in MHD disk 
turbulence have a typical coherence time of the order of the orbital period, and as a 
consequence are able to exchange angular momentum with an embedded planet across a 
range of disk radii (not only at narrow resonances). The study by Nelson & Papaloizou 
(2004) was based on both global ideal MHD disk simulations, with an aspect ratio of 
h/r = 0.07, and local shearing box calculations. The global runs realized an effective 
Shakura-Sunyaev a = 7 x 10~^, which if replicated in a real disk would be consistent 
with observational measures of T Tauri disk evolution (Hartmann et al. 1998). For all 
masses considered in the range 3 ^ Mp 30 M^, the instantaneous torque on the 
planet from the MHD turbulent disk exhibited large fluctuations in both magnitude and 
sign. Averaging over ss 20 orbital periods, the mean torque showed signs of converging 
to the Type I rate, although the rate of convergence was slow, especially for the low- 
est mass planets in the global runs. These properties are generally in accord with other 
studies of the variability of MHD disk turbulence (Hawley 2001; Winters, Balbus & Haw- 
ley 2003a). Very roughly, the Nelson & Papaloizou (2004) simulations suggest that up 
to Mp ~ 10 the random walk component dominates steady Type I drift over time 
scales that substantially exceed the orbital period. 

We caution that existing studies of the stochastic migration regime are unrealistic. 
Ideal MHD is not a good approximation for the protoplanetary disk at those radii where 
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planet formation occurs, and there may be dead zones in which MHD turbulence and 
angular momentum transport is highly suppressed (Gammie 1996; Sano et al. 2000; Fro- 
mang, Terquem & Balbus 2002; Salmeron & Wardlc 2005). Wc also note that a significant 
random migration component, if it indeed adds to rather than supplanting steady Type I 
migration, does nothing (on average) to help the survival prospects of low mass planets. 
Nevertheless, if turbulent fluctuations (whatever their origin) do occur in the disk, the 
resulting random walk migration could be important for planet formation. In the terres- 
trial planet region, stochastic migration might deplete low mass planetary embryos that 
would be relatively immune to ordinary Type I migration, while simultaneously promot- 
ing radial mixing and collision of planetesimals. For giant planet formation, a significant 
random component to core migration would have the effect of creating large fluctuations 
in the planetesimal accretion rate, while also potentially acting to diffuse the planetesi- 
mal surface density. Large fluctuations in the planetesimal accretion rate favor the early 
onset of rapid gas accretion, and allow for the flnal core mass to be substantially smaller 
than would be expected in the case of a static core (Rice & Armitage 2003). 



4. Type II migration 

4.1. Conditions for the onset of Type II migration 

In a viscous disk, the threshold between Type I and Type II migration can be derived by 
equating the time scale for Type I torques to open a gap (in the absence of viscosity) with 
the time scale for viscous diffusion to fill it in (Goldreich & Tremaine 1980; Papaloizou & 
Lin 1984). For a gap of width Ar around a planet with mass ratio q = Mp/M^, orbiting 
at distance Vp, Type I torques can open the gap on a timescale (Takeuchi, Miyama & 
Lin 1996) 

"^^P^" ^ w?q^ \rp) P • 
Viscous diffusion will close the gap on a time scale, 

'^close ) 
U 

where v, the kinematic viscosity,is usually written as v = aCg/fij, (Shakura & Sunyaev 
1973). Equating Topcn and Tdosc, and noting that waves with m w rp/h dominate the 
Type I torque, the condition for gap opening becomes, 

.a(i)>-. 

For h/r ~ 0.05 and a = 10~^, the transition (which simulations show is not very sharp) 
occurs at Qcrit ~ 2.5 x 10~^, i.e. close to a Saturn mass for a Solar mass star. Since the 
most rapid Type I migration occurs when q w gcrit, we can also estimate a minimum 
migration time scale by combining the above expression with the time scale formula of 
Tanaka, Takeuchi & Ward (2002). This yields, 

r^i„~(2.7+l.l/3)-i||a-V2fi-i, 



and is almost independent of disk properties other than the local mass. 
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4.2. The rate of Type II migration 

Once a planet has becomes massive enough to open a gap, orbital evolution is predicted 
to occur on the same local time scale as the protoplanetary disk. The radial velocity of 
gas in the disk is, 

M 



27rrS ' 

which for a steady disk away from the boundaries can be written as, 

2 r 

If the planet enforces a rigid tidal barrier at the outer edge of the gap, then evolution of 
the disk will force the orbit to shrink at a rate rp '^Vr, provided that the local disk mass 
exceeds the planet mass, i.e. Trr^E > Mp. This implies a nominal Type II migration time 
scale, valid for disk dominated migration only. 

For h/r = 0.05 and a = 10~^, the migration time scale at 5 AU is of the order of 0.5 Myr. 

In practice, the assumption that the local disk mass exceeds that of the planet often 
fails. For example, a /3 = 1 disk with a mass of 0.01 Mq within 30 AU has a surface 
density profile, 



2 



The condition that TrrpS = Mp gives an estimate of the radius within which disk domi- 
nation ceases of. 

Interior to this radius, the planet acts as a slowly moving barrier which impedes the 
inflow of disk gas. Gas piles up behind the barrier - increasing the torque - but this 
process does not continue without limit because the interaction also deposits angular 
momentum into the disk, causing it to expand (Pringle 1991). The end result is to slow 
migration compared to the nominal rate quoted above. For a disk in which the surface 
density can be written as a power-law in accretion rate and radius, 

S oc MV, 

Syer & Clarke (1995) define a measure of the degree of disk dominance. 



B = 



— p 



Mp 



Then for B < 1 (the planet dominated case appropriate to small radii) the actual Type II 
migration rate is (Syer & Clarke 1995), 

Note that with this definition of B, disk dominance extends inward a factor of a few 
further than would be predicted based on the simple estimate given above. 

Evaluating how the surface density depends upon the accretion rate - and thereby 
determining the a which enters into the suppression term - requires a full model for the 
protoplanetary disk (not just knowledge of the slope of the steady-state surface density 
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profile). For the disk models of Bell et al. (1997), we find that a ~ 0.5 at 1 AU for 
M ~ 10~^ M0yr~-^. At this radius the model with a = 10~^ has a surface density of 
about 200 gcm~^. For a Jupiter mass planet we then have B = 0.3 and t// = I.Stq. This 
is a modest suppression of the Type II rate, but the effect becomes larger at smaller radii 
(or for more massive planets). It slows the inward drift of massive planets, and allows a 
greater chance for them to be stranded at sub-AU scales as the gas disk is dissipated. 

These estimates of the Type II migration velocity asstune that once a gap has been 
opened, the planet maintains an impermeable tidal barrier to gas inflow. In fact, simu- 
lations show that planets are able to accrete gas via tidal streams that bridge the gap 
(Lubow, Seibert & Artymowicz 1999). The effect is particularly pronounced for planets 
only just massive enough to open a gap in the first place. If the accreted gas does not 
have the same specific angular momentum as the planet, this constitutes an additional 
accretion torque in addition to the resonant torque. It is likely to reduce the Type II 
migration rate further. 

4.3. Simulations 

Simulations of gap opening and Type II migration have been presented by a large num- 
ber of authors, with recent examples including Bryden et al. (1999), Lubow, Seibert & 
Artymowicz (1999), Nelson et al. (2000), Klcy D'Angelo & Henning (2001), Papaloizou, 
Nelson & Masset (2001), D'Angelo, Henning & Kley (2002), D'Angelo, Kley & Henning 
(2003), D'Angelo, Henning & Klcy (2003), Bate et al. (2003), Schafor et al. (2004) and 
Lufkin et al. (2004). These authors all assumed, for simplicity, that angular momentum 
transport in the protoplanetary disk could be represented using a microscopic viscosity. 
Only a few recent simulations, by Winters, Balbus & Hawley (2003b), Nelson & Pa- 
paloizou (2004) and Papaloizou, Nelson & Snellgrove (2004), have directly simulated the 
interaction of a planet with a turbulent disk. For planet masses significantly above the 
gap opening threshold, simulations support the general scenario outlined above, while 
also finding: 

1. Significant mass accretion across the gap. For planet masses close to the gap 
opening threshold, accretion is surprisingly efficient, with tidal streams delivering gas at 
a rate that is comparable to the disk accretion rate in the absence of a planet (Lubow, 
Seibert & Artymowicz 1999). It is also observed that accretion cuts off rapidly as the 
planet mass increases toward 10 Mj, giving additional physical motivation to the stan- 
dard dividing line between massive planets and brown dwarfs. 

2. Damping of eccentricity. Goldreich & Tremaine (1980) noted that the evolution 
of the eccentricity of a planet embedded within a disk depends upon a balance between 
Lindblad torques, which tend to excite eccentricity, and corotation torques, which damp 
it. Recent analytic work (Ogilvie & Lubow 2003; Goldreich & Sari 2003; and references 
therein) has emphasized that if the corotation resonances arc even partially saturated, 
the overall balance tips to eccentricity excitation. To date, however, numerical simula- 
tions (e.g. those by Papaloizou, Nelson & Masset 2001) exhibit damping for planetary 
mass bodies, while ciccentricity growth is obtained for masses Mp > 20 Mj in the brown 
dwarf regime. In a recent numerical study, however, Masset & Ogilvie (2004) present ev- 
idence that the resolution attained by Papaloizou, Nelson & Masset (2001) was probably 
inadequate to determine the sign of eccentricity evolution for Jupiter mass planets, which 
remains uncertain. Further numerical work is needed. 

4.4. Comparison with statistics of extrasolar planetary systems 

The estimated time scale for migration of a giant planet from 5 AU to the hot Jupiter 
region is of the order of a Myr. This time scale is short enough - compared to the lifetime 
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semi— major axis / AU a / AU 

Figure 4. Left hand panel: the distribution of extrasolar planets discovered via radial velocity 

surveys in the a-Mp sin(i) piano. The dashed diagonal traces a line of equal case of dctcctability - 
planets on circular orbits lying along lines parallel to this cause the same amplitude of stellar ra- 
dial velocity variations. Right hand panel: the number of (lct(>(:t(xi ])lancts with Alp sin(i) > 1 A4j 
is plotted as a function of semi-major a^xis. The solid curve shows the predicted distribution from 
a pure migration model by Armitage et al. (2002). The theoretical curve is unaltered from the 
2002 version except for an arbitrary normalization. 



of typical protoplanetary disks - to make migration a plausible origin for hot Jupitcrs, 
while not being so short as to make large-scale migration inevitable (the latter would 
raise obvious concerns as to why there is no evidence for substantial migration of Jupiter 
itself). Having passed this initial test, it is then of interest to try and compare quantitative 
predictions of migration, most obviously the expected distributions of planets in mass 
and orbital radius, with observations. Models that attempt this include those by Trilling 
et al. (1998), Armitage et al. (2002), Trilling, Lunine fc Benz (2002), Ida & Lin (2004a) 
and Ida & Lin (2004b) . Accurate knowledge of the biases and selection function of the 
radial velocity surveys is essential if such exercises are to be meaningful, making analyses 
such as those of Gumming, Marcy & Butler (1999) and Marcy et al. (2005) extremely 
valuable. 

Figure 4 shows how the distribution of known extrasolar planets with semi-major axis 
compares to the pure migration model of Armitage et al. (2002). In this model, we 
assumed that giant planets form and gain most of their mass at an orbital radius (specif- 
ically 5 AU) beyond where most extrasolar planets are currently being detected. Once 
formed, planets migrate inward via Type II migration and are either (a) swallowed by the 
star, or (b) left stranded at some intermediate radius by the dispersal of the protoplane- 
tary disk. We assumed that disk dispersal happens as a consequence of photoevaporation 
(Johnstone, Hollenbach & Bally 1998), and that the probability of planet formation per 
unit time is constant over the (relatively short) window during which a massive planet 
can form at 5 AU and survive without being consumed by the star. Although clearly 
oversimplified, it is interesting that this model continues to reproduce the orbital dis- 
tribution of known planets out to radii of a few AU, once selection effects have been 
taken into account (here, by simply ignoring low mass planets that are detectable only 
at small orbital radii). Moreover, it predicts that substantial outward migration ought to 
occur in disks where strong mass loss prompts an outwardly directed radial velocity in 
the giant planet forming region (Veras & Armitage 2004). Planets at these large radii are 
potentially detectable today via their effect on debris disks (Kuchner & Holman 2003). 
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Figure 5. The minimum metallicity required to form a gas giant planet as a function of orbital 
radius, from simplified core accretion calculations by Rice & Armitage (2005). The host metal- 
licity is expressed via the parameter /dust, which is proportional to lO''^'^^^'. The models assume 
that dispersion in disk metallicity dominates over dispersion in either disk gas mass or disk 
lifetime in controlling the probability of planet formation. Type I migration of the core prior 
to accretion of the envelope is included, using several different prescriptions indicated by the 
different symbols, though none of the Type I rates are as large as the baseline rate of Tanaka, 
Takeuchi & Ward (2002). 



Additional clues to the role of migration in forming the observed population of extra- 
solar planets may be possible by combining largo planet samples with knowledge of the 
host stars' metallicity. It is now clear that the frequency of detected planets increases 
rapidly with host metallicity, and that the measured metallicity reflects primarily the 
primordial composition of the star-forming gas rather than subsequent pollution of the 
convective envelope (Santos, Israelian & Mayor 2004; Fischer &; Valenti 2005; and refer- 
ences therein). The existence of such a correlation is not in itself surprising, given that the 
time scale for core accretion decreases quickly with increasing surface density of planetes- 
imals. However, the sharpness of the rise in planet frequency over a fairly narrow range of 
stellar [Fe/H] is striking, since it suggests that metallicity, rather than variations in ini- 
tial gas disk mass or gas disk lifetime, may well be the single most important parameter 
determining the probability of giant planet formation around a particular star. 

Motivated by these observations, we have investigated the use of the critical or thresh- 
old metallicity for giant planet formation as a discriminant of different planet formation 
models (Rice & Armitage 2005). Using simplified core accretion models similar to those 
of Ida & Lin (2004a), we have calculated the radial dependence of the threshold metal- 
licity under the assumption that disks around different stars have similar initial gas 
masses and lifetimes (this could be replaced with the much weaker assumption that the 
gas mass and disk lifetime are not correlated with the metallicity). A sample result is 
shown in Figure 5. By definition, planets that just manage to form as the gas disk is 
being lost suffer little or no Type II migration, so delineating the threshold metallicity 
curve observationally can yield information on planet formation that is independent of 
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Type II migration uncertainties. We find that the most important influence on the shape 
of the threshold metallicity curve is probably Type I core migration prior to accretion 
of the gas envelope. When this is included we derive a monotonically rising minimum 
metallicity curve beyond about 2 AU. In the absence of significant core migration, the 
threshold metallicity is flat beyond the snow line (with a weak dependence on the surface 
density profile of planetesimals), and the location of the snow line; may be preserved in 
the observed distribution of planets in the orbital radius / metallicity plane. 

5. Summary 

Migration via gravitational interactions between planets and the gaseous protoplane- 
tary disk appears to be central to understanding both the formation of gas giant planets, 
and their early orbital evolution to yield the cxtrasolar planetary systems currently being 
observed. Although there are uncertainties in our understanding of migration, there has 
also been enough progress to convince us that Type I migration is probably a vital ingre- 
dient in the formation of gas giant planets via core accretion. This is trivially true using 
the current best estimates of the Type I migration time scale, but it would still remain 
important even if the rate was suppressed by as much as two orders of magnitude. On a 
similarly firm footing is the assertion that gas disk migration - probably in the Type II 
regime - is responsible for the existence of most of the hot Jupiters. Although other mi- 
gration mechanisms certainly exist, it requires moderate care to avoid substantial orbital 
evolution once a planet has formed in a gas disk. 

Equally interesting are the major unknowns. Does turbulence within the disk lead to 
random walk migration of low mass bodies, and if so, is this important for terrestrial 
planet formation and / or core accretion? Do corotation torques qualitatively change the 
behavior of migrating planets with masses just above the gap-opening threshold? Is the 
eccentricity of massive planets excited by the interaction with the gas disk? A positive 
answer to any of these questions would require substantial changes to our ovca'all picture 
of planet formation. Addressing them will probably require, in part, high resolution 
simulations that include more of the complex physics of angular momentum transport 
within the protoplanetary disk. 
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